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We present a systematic study of the photo-absorption spectra of various Si„Hm clusters (n=l-10, m=l- 
14) using the time-dependent density functional theory (TDDFT). The method uses a real-time, real-space 
implementation of TDDFT involving full propagation of the time dependent Kohn-Sham equations. Our results 
for SiH4 and Si2Hg show good agreement with the earlier calculations and experimental data. We find that for 
small clusters (n<7) the photo-absorption spectrum is atomic-like while for the larger clusters it shows bulk-like 
behaviour. We study the photo-absorption spectra of silicon clusters as a function of hydrogenation. For single 
hydrogenation, we find that in general, the absorption optical gap decreases and as the number of silicon atoms 
increase the effect of a single hydrogen atom on the optical gap diminishes. For further hydrogenation the 
optical gap increases and for the fully hydrogenated clusters the optical gap is larger compared to corresponding 
pure silicon clusters. 
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I. INTRODUCTION 

Recently, there has been a renewed interest in understand- 
ing the optical properties of clusters because confinement of 
electrons changes the physical properties. Thus by varying 
size of the clusters the optical properties can be tuned accord- 
ing to the desired application^ In particular, optical proper- 
ties of silicon and hydrogenated silicon clusters have been of 
great interest due to the observation of photo luminescence 
(PL) in porous silicon"'^. The structure and properties of sili- 
con clusters can be tuned by varying the cluster size as well as 
doping. An important dopant for silicon clusters is hydrogen 
and it plays an important role in structural stability. Experi- 
mental studies also confirm this fact"*^, but in spite of several 
investigations many issues about hydrogenated silicon clus- 
ters have not been understood. It is not clear how the structure 
and optical properties of the cluster evolves with size and as 
a function of hydrogenation. It is therefore interesting to cal- 
culate the photo-absorption (R\) spectra and compare it with 
experiment. 

There has been a considerable amount of work experimen- 
tal as well as theoretical on silicon and hydrogenated silicon 
clusters. Suto and Lee measured the photo-absorption and flu- 
orescence of silane in the energy range 8-12 eV using syn- 
chrotron radiation^. The optical absorption of silane and dis- 
ilane has been measured by Itoh et. al. in the energy range 
6-12 eV-. Cheshnovsky et. al. have measured the photo- 
electron spectra of charged silicon and germanium clusters^. 
Rinnen and Mandich have measured the photo-dissociation 
spectra of neutral silicon clusters Si„(n=18-41)^. Murakami 
and Kanayama investigated the stability of some silicon and 
hydrogenated silicon clusters using a quadruple ion trapi. 
More recently Antonietti et. al. deduced the photo-absorption 
spectra of charged silicon clusters from photo-dissociation 
of charged xenon-silicon clusters^. On the theoretical side 
Chantranupong et. al. performed a configuration interac- 
tion (CI) calculation for a large number of low lying states 



in silane^". Rubio et. al. calculated the photo-absorption 
spectra of silicon and alkali metal clusters using time de- 
pendent local density approximation (TDLDA)". Rohlfing 
and Louie calculated the optical absorption spectra of hydro- 
gen terminated silicon clusters by solving the Bethe-Salpeter 
equationi^. Vasiliev et. al.'^ and Marques et. al. ^^calculated 
the optical absorption spectra of Si„H„, clusters using lin- 
ear response theory within TDLDA. Rao et. al. measured 
the photo luminescence of a dispersion of 1 nm silicon par- 
ticles obtained from crystalline silicon that is dispersed into 
nanoparticles through electrochemical etching with HF and 
H202'^. They also calculated the photo-absorption spectra 
of Si29H24 using time dependent density functional theory 
(TDDFT). Lehtonen and Sundholm calculated the absorption 
spectra of three hydrogen terminated silicon clusters using 
TDDFTi^. 

The earlier studies on the optical properties of silicon and 
hydrogenated silicon clusters focused on the size dependence 
of the PL and photo-absorption^si^"^. Some of these studies 
ignored the influence of oscillator strength's (electric dipole 
matrix elements) and hence were not in good agreement with 
the experimental data. Calculations based on DFT using LDA 
and the generalised gradient approximation (GGA) which in- 
cluded the dipole matrix elements suffered from the drawback 
that is inherent in LDA/GGA i.e. the energy gaps were un- 
derestimated because these calculations ignored the effect of 
excited states. To overcome this drawback, the configuration 
interaction method or the methods based on solving the Bethe- 
Salpeter equation along with the GW approximation have 
been suggested'-. These methods require a lot of computer 
time and hence have been restricted to small clusters. A com- 
putational technique based on linear response theory within 
TDLDA had been proposed by Vasiliev et. al.'^. This is a 
natural extension of the LDA ground state density functional 
formulation designed to include excited states. This method 
is faster than the BS and GW methods and can therefore be 
used for larger clusters. Vasiliev et. al. have shown the vi- 



ability of their method by performing calculations on silicon 
and hydrogenated and oxygenated silicon clusters. Another 
implementation of TDDFT has been formulated by Castro et. 
al. which allows the calculation of the excited state energies 
and optical absorption spectra. In the present work, we have 
used this method^. 

Most of earlier calculations have considered only hydrogen 
terminated silicon clusters. There seems to be a lack of a sys- 
tematic study on small hydrogenated silicon clusters. In this 
paper we report calculations for the silicon and hydrogenated 
silicon clusters in addition to some hydrogen terminated sili- 
con clusters. Thus our calculations will show the effect of hy- 
drogenation on the optical properties of silicon clusters. Our 
emphasis is on the smaller clusters so as to bring out the evo- 
lution of the optical properties as we increase the number of 
silicon and hydrogen atoms in the cluster 

The plan of the paper is as follows. In section II, we briefly 
discuss the method and give computational details. In section 
III we present our results and discussion. In section IV we 
present our conclusions. 



II. METHOD AND COMPUTATIONAL DETAILS 

We have used TDDFT for our calculation of the photo- 
absorption spectrum^^'^"*^. For the sake of completeness we 
shall summarise the essentials of the method. In TDDFT, the 
basic variable is the one electron density n(r,f), which is ob- 
tained with the help of a fictitious system of non-interacting 
electrons, the Kohn-Sham system. The time-dependent Kohn- 
Sham equations are 



function 5(a)) 



/— V/;(r,0 - [-VV2 + v^5(r,f)] v/,(r,0 



dt 



(1) 



where \j/i{r,t) are Kohn-Sham one electron orbitals. In terms 
these orbitals n(r,f ) can be written as 



«(r,f) 
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(2) 



The Kohn-Sham potential can be written as 

VKs{r,t) = Va,{r,t) + VHartreeir,t) +V_^cir,t) (3) 

where the first term is the external potential, the second 
Hartree potential and the last the exchange and correlation 
potential. For obtaining this potential we use adiabatic local 
density approximation (ALDA). 

In our work we shall use the TDDFT scheme which cal- 
culates propagation of the time-dependent Kohn-Sham equa- 
tions in real time. In this scheme the electrons are given some 
small momentumCfc) to excite all the frequencies^^. This is 
achieved by transforming the ground state wave function ac- 
cording to 



VA,(r,5f)=e''^->,(r,0) 



(4) 



and then propagating these wave functions for some finite 
time. The spectrum is then obtained from dipole strength 



S{(o) = — Ima{co) 

K 



(5) 



where the a(a))is the dynamical polarizability and is given by 



a{(o) = - I dte'""[d{t)-d{0)] 



(6) 



where d{t) is dipole moment of the system. 

We can also define a quantity know as the oscillator strength 
to express the strength of the transition. 
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where 0o and 0/ are the ground and excited state respectively. 
The oscillator strength is related to the dipole strength func- 
tion defined earlier using the following relationship 



S{(0)=Y,f[Xd{(0-COi) 

I 



(8) 



The sum over the oscillator strength gives the number off ac- 
tive electrons in the system. 



Lf'-N 



(9) 



where A^ is the number of active electrons in the system. 

The initial structures of the clusters used for our calcula- 
tions have been obtained by Balamurugan and Prasad^^^ ear- 
lier from their analysis based on the Car-Parrinello molecular 
dynamics (CPMD)^**. The resulting structures have been fur- 
ther optimised using the electronic structure method imple- 
mented in VASP (Vienna Ab-initio Simulation Package)^'^ . We 
do not find any significant changes in the structures compared 
to the CPMD results. In fact the starting CPMD configura- 
tion gives forces that are around 0.5 eV/A which eventually 
reduce to around 0.05 eV/A in VASP. VASP employs den- 
sity functional theory (DFT) and we have employed the lo- 
cal density approximation (LDA) for the exchange correlation 
using norm conserving pseudo-potentials. The optimisation is 
done only by relaxing the ions via the conjugate gradient (CG) 
method and using a k-point Monkhorst-Pack mesh of 4x4x4. 
All calculations have been performed in a cubic supercell of 
length 20 A. 

In this work the photo-absorption spectra of the optimised 
structures was studied using OCTOPUS code^^-^" where the 
above approach is implemented. All calculations are ex- 
panded in a regular mesh in real space, and the simulations 
are performed in real time. The local density approximation 
was employed to keep consistency with the geometry optimi- 
sation process. OCTOPUS uses a uniform grid in real space, 
which is located inside the sum of n spheres, one around each 
atom of the n-atom cluster. For all clusters a minimisation of 
energy with respect to the radius and the grid spacing was car- 
ried out. We required the radius of each sphere to be 6-8 A and 
the grid spacing 0.28-0.4 A for optimal energy minimisation. 

For the TDDFT the propagation in real time has been per- 
formed with 30,000 time steps with a total simulation time 



of around 124 fs. This gives a good resolution of the spec- 
tra. Throughout the calculations the ions were kept static and 
in order to approximate the evolution operator the approxi- 
mated enforced time reversal symmetry (aetrs) method was 
employed. Numerically the exponential of the Hamiltonian, 
which is used to approximate the evolution operator was eval- 
uated using a simple Taylor expansion of the exponential. 



III. RESULTS AND DISCUSSIONS 

In order to verify our calculations with the experimental 
data, we chose to compare the two most stable clusters SiH4 
and Si2H6 with the available experimental data as our bench- 
mark. In optical properties we are interested in transitions 
form the occupied levels to the unoccupied levels. The struc- 
tures in the photo-absorption spectra are identified with these 
transitions. It is therefore interesting to calculate these energy 
differences and compare them with the energy differences de- 
duced from the experimental spectra. Table I shows such a 
comparison. For our calculations, we have checked that the 
f-sum rule is satisfied^' . In Table 1 we have also included the 
Bethe-Salpeter results of Rohlfing and Louie'^, the TDLDA 
results of Vasilev et. alji^ and the TDLDA results from Oc- 
topus of Marques et. al4^ along with the transitions identifi- 
cation. The uppermost occupied states result from a hybridi- 
sation of the silicon and hydrogen states while the lowest un- 
occupied states are primarily silicon states. Thus the energy 
gap is dependent on the the bonding and anti-bonding silicon 
states. The transition between these states have been identi- 
fied as 4s, 4p, 4d states (these refer to the angular momentum 
character of the final states). We find good agreement with 
the experimental data of Itoh et. al.^and Suto and Lee^ and 
other theoretical calculations. In particular, we get very good 
agreement with the results of Marques et. al.''*. As compared 
to the work of Marques et. al. we obtain an additional peak 
at 8.3 eV for Si2H6 which is observed in other worksi^ and 
the experiments^'^. Our simulations pick up this additional 
peak because our broadening factor, which is the inverse of 
the time steps in the TDDFT run is quite narrow as compared 
to Marques et. al. Fig. 1 shows the structures and photo- 
absorption spectra of the current work, Vasiliev et.al.^^ and 
Marques et. al.''* for SiH4 and Si2H6 clusters. The calculated 
photo-absorption cross-section (using TDDFT) seems to be in 
good agreement with the TDLDA results of Vasiliev et. alji2, 
and the TDLDA results of Marques et. al.''* using Octopus. 

In order to see the effect of increasing the size of silicon 
clusters on the optical spectra, we present in Fig. 2 I the opti- 
mised structures of silicon clusters Si„ [n=l-10] along with the 
photo-absorption spectra. To see the effect of adding a single 
H atom to these clusters, we have also calculated the photo- 
absorption spectra of Si„H clusters which are also shown in 
Fig. 2 II along with the silicon clusters results. Each figure 
shows the structure of silicon and hydrogenated silicon clus- 
ters and the calculated photo-absorption cross-section. Con- 
sider first the silicon clusters. For small clusters (up to n=7) 
we find that the photo-absorption spectra is a combination of 
many peaks and looks like that of isolated atoms. However 



Figure 1: Structure and photo-absorption spectra of present work 
(Black solid line), Vasiliev et.al. (Blue dashed line) and Marques 
et. alili (Red solid line) for SiH4 and Si2H6 clusters. 
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Figure 2: I: Structure and photo-absorption spectra of Si„ clusters 
(n=l-10 from top to bottom) ; II: Structure and photo-absorption 
spectra of Si„H clusters (n=l-10 from top to bottom). 
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for n>7, the optical spectra looks bulk-like^^. This can be 
understood using a simple tight-binding picture. In a larger 
silicon cluster, the overlap between electronic wave-functions 
lifts the degeneracy of the energy levels resulting in bunch- 
ing of energy levels in a narrow energy range. This results 
in broadening of PA spectrum for larger clusters. In Fig 2 
we see that the main structure in the PA spectra is located at 
around 9 eV and a minor structure starts to build up at around 
15 eV. The same trends are found in the singly hydrogenated 
clusters. For smaller clusters (n<7) the PA spectrum changes 
significantly as we increase the number of hydrogen atoms, 
while for larger clusters this change is small. 

To see the effect of further hydrogenation, we present the 
PA spectra of hydrogenated silicon clusters where the number 
of hydrogen atoms is increased in Fig. 3. We have performed 



Table I: Excitation energies in eV for SiH4and Si2H6 clusters. 
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calculations for SisHs, Si4H4 and Si4H8, SisHg and Si5Hi2 
and SieH? and SieH^ clusters. The different positions of the 
hydrogen atoms in SisH and Si3H3 makes the two structures 
very different. Hence there is a large difference between the 
PA spectra. In going from SisH to Si3H3 we find that the 
number of structures in the PA spectra increases and some of 
the new structures have larger peaks. In fact the PA spectrum 
goes from atom-like to bulk-hke. We see similar trends in 
Si4H4 and Si4H8 except that the siUcon backbone is a rhom- 
bus. We find that the new structures in Si4H8 PA spectrum 
have significantly larger peaks compared to Si4H4. In these 
hydrogenated clusters (SisHg, 814114, Si4H8, SisHg , Si5Hi2, 
SieHv and Si6Hi4) which have different silicon backbones, 
we observe that the bulk-like behaviour seems to start at n=4 
which is much earlier than in singly hydrogenated and silicon 
clusters. This effect may be attributed to the large number of 
hydrogen atoms. 

We obtained two very closely related structures of SijH 
having a difference of about 0.043 eV in their total energy. 
Both the structures have been shown in Fig. 4 and the cor- 
responding PA spectra in Figure 2 II. The silicon backbone in 
both these clusters remains the same and the only difference is 
found in the position of the hydrogen atom. In the ground state 
structure (solid line in the PA spectra) hydrogen is bonded to 
two silicon atoms while in the other (dashed line in the PA 
spectra) it is bonded to only one silicon atom. The PA spectra 
of the two structures is quite different. This shows that the PA 
spectrum is sensitive to small changes in the structure. Thus 
the PA spectrum can be used to identify the structure of the 
clusters which differ by a very small energy. 

In Table II we present the optical absorption gaps of the var- 
ious clusters. We define the optical gaps through the integral 
oscillator strength rather than as the energy of the first dipole 
allowed transition in the absorption spectra. The integral os- 
cillator strength gives the total number of active electrons in 
the system. In this approach the value of the optical absorp- 
tion gap is determined at a very small but non-zero fraction 
of the complete oscillator strength^^. We set this threshold to 
10^"* of the total oscillator strength. This value is chosen be- 
cause it stands above the value of "numerical noise"and at the 
same time it is sufficiently small so as to not suppress the ex- 
perimentally detectable dipole allowed transitions. This def- 
inition of absorption gap doesn't affect the values of optical 
gaps for small clusters, since the intensity of first transitions 
is well above the selected threshold. 



Figure 3: Structure and photo-absorption spectra of Si„H„, clusters 
(from top to bottom: Si3H3, Si4H4, Si4H8, SisHj,, Si5Hi2, SigHy, 
Si6Hi4) 
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Figure 4: Structures of ground state (I) and higher energy state (II) 
SisH clusters. 
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In Fig. 5 we plot the optical gap as a function of the num- 
ber of silicon atoms (n) for various Si„ (solid Line) and Si„H 
(dashed Line) clusters (n=l-10). It shows a general trend that 
the addition of a single hydrogen atom reduces the optical gap 
as compared to the silicon clusters. As the number of silicon 
atoms in the cluster increase we find that the effect of a sin- 
gle hydrogen atom reduces the optical gap and the difference 
between the optical gap of Si„ and Si„H gradually decreases 
with the increase in silicon atoms. In the bulk limit a single 
hydrogen atom should not distort the optical gap and thus in 
that limit the optical gap of the silicon will be the same as the 
optical gap of the singly hydrogenated silicon. 

To see the effect of further hydrogenation, we list the op- 



Figure 5: Optical absorption gap of various Si„ (solid line filled cir- 
cle) and Si„H (dashed line empty circle) clusters as a function of the 
number of silicon atoms (n). 



Table II: Optical gaps of various clusters. 
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tical gaps of all the clusters in Table II. We see that addition 
of a single hydrogen atom reduces the optical gap, in general, 
but with further hydrogenation the optical gap increases. For 
the fully hydrogenated clusters the optical gap is larger than 
the unhydrogenated clusters. Thus in the bulk limit the hy- 
drogenated system will show a larger gap which is consistent 
with the experimental observationsi^ 



IV. CONCLUSIONS 

The photo-absorption spectra of silicon and hydrogenated 
silicon clusters have been calculated using TDDFT. Our cal- 
culations show very good agreement for the benchmarks of 
SiH4 and Si2H6 with the earlier theoretical calculations and 
experimental results. We find that the changes in the photo- 
absorption spectra correlates well with the changes in the 
structure of the cluster. In the singly hydrogenated clusters 
we find bulk-like behaviour for larger clusters (n>7) while for 
the smaller clusters we find the PA spectra is composed of 
numerous peaks as in atoms. We find that the PA spectra of 
two structures of SijH which are close in energy to be quite 
different. For the other hydrogenated clusters the bulk-like 
behaviour appears for smaller n which can be attributed to 
the larger number of hydrogen atoms. The addition of a sin- 
gle hydrogen to the silicon cluster shows a decrease in the 
absorption optical gap. The optical gap decreases gradually 
with the number of silicon atoms in the cluster giving the bulk 
limit in which an addition of a single hydrogen to the silicon 
should not effect the optical gap of silicon. With further hy- 
drogenation the optical gap increases and for the fully clusters 
the optical gap is larger than the corresponding gap for pure 
silicon clusters. 
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